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ABSTRACT 

We present Hubble Space Telescope near-infrared transit photometry of the nearby hot- Jupiter HD189733b. The obser- 
vations were taken with the NICMOS instrument during five transits, with three transits executed with a narrowband 
filter at 1.87 /im and two performed with a narrowband filter at 1.66 /im. Our observing strategy using narrowband 
filters is insensitive to the usual HST intra-orbit and orbit-to-orbit measurement of systematic errors, allowing us to 
accurately and robustly measure the near-IR wavelength dependance of the planetary radius. Our measurements fail to 
reproduce the Swain et al. absorption signature of atmospheric H2O below 2(im at a 5a confidence level. We measure 
a planet-to-star radius contrast of 0.15498T0. 00035 at 1.66 /im and a contrast of 0.15517±0. 00019 at 1.87 fj,m. Both 
of our near-IR planetary radii values are in excellent agreement with the levels expected from Rayleigh scattering by 
sub-micron haze particles, observed at optical wavelengths, indicating that upper-atmospheric haze still dominates the 
near-IR transmission spectra over the absorption from gaseous molecular species at least below 2 \im. 
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1. Introduction 

Transiting close-in exoplanets have provided the extraordi- 
nary possibility to begin studying the detailed character- 
istics of extrasolar planets. A transit or anti-transit event 
allows for an exoplanet to be temporally resolved from the 
bright parent host star, isolating such valuable information 
as the planet's absorption or emission spectra. During a 
primary transit event, both the opaque body of the planet 
as well as its atmosphere blocks light from the parent star. 
A precise determination of the radius of the planet can be 
made from the total obscuration, while partial transmission 
of light through the exoplanet's atmosphere, with its wave- 
length and altitude dependence, allows for detecting com- 
position and structure. The frequent transits, anti-transits, 
and large signals make hot Jupiters the best targets for 
these studies. Space-based observatories, notably HST and 
Spitzer, have proven to be extremely efficient platforms on 
which to do follow-up transit and anti-transit studies given 
their superior photometric performance, wide wavelength 
range, and ability to observe multiple transits over long 
baselines within a single observing season. 

Two hot Jupiters in particular, HD189733b and 
HD209458b, currently offer the very best laboratories in 
which to study exoplanet atmospheres and have become the 
prototype hot- Jupiter planets. These two planets have the 
brightest parent stars among transiting planets and have 
large transit depths, making precise studies at high signal- 



to-noise ratios possible. The first transiting planet discov- 
ered, HD209458b, holds the distinction of the first detection 
of an extrasolar planetary atmos phere (ICharbonneau et al] 
2002) and escaping atmosphere (IVidal-Madiar et aLl l2003, 
12004 120081 ; iBallester etall |2007j). The optical transmis- 
sion spectra of this planet shows ev idence for several dif- 
ferent layers of Na (jSing et al.ll2008aH bl) as well as Rayleigh 
scatterin g by molecular hydrogen (X ecavelie r des Etangs 
et_al. [2008bJ) and the presence of TiO/VO (jDesert et al.1 
l2008f ). The atmospheric Na signature has also been con - 
firmed by ground based observations (|Snellen et aLl l^OOS). 
This planet also exhibits a stratospheric temper ature inver- 
sion (|Burrows et al.ll2007l ; iKnutson et aTll2008l ) thought to 
be c aused by a strong optical absorber such as TiO f Desert 
et al. l2008tlHubenv et al.l|2003l ; iFortnev et al.ll2008T l and/or 
W 



HS and S 2 (jZahnle et al. 2009). 
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HD189733b (Bou chv et aT]|2005l: Hebrard & Lecavelier 
Des Etangs |2006|) is among the closest known transiting 
planets with a K1.5V type parent star, giving it one of the 
largest transit and anti-transit signals known. Spitzer anti- 
transit measurements have shown efficient heat redistribu- 
tion, measuring th e planet's temperature pr ofile from or- 
bital phase curves (|Knutson et al.ll2007ll2009f ) and a defini- 
tive detection of atm ospheric water from emission spectra 
(|Grillmair et al.|[2008l ). The signature of CO2 has also been 
seen from HST/NICMOS emission spectra (|Swain et al.l 
120091 ). In primary transit. IPont et all (|2008l ) used the ACS 
grism to provide the first transmission spectra from 0.6 to 
1.0 /Ltm. This spectrum is seen to be almost featureless likely 
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indicating the presence of high altitude haze, with a A~ 4 
wavelength dependence of the spectra likely due to Raylcigh 
scattering b y sub- micro n MgSi Os molecules (Lecavelier des 
Etangs et al.|2008a|). From the ground, |RedfielcUtaL| (|2008l ) 
detected stro ng Na absorp t ion in the core of the doublet. In 
the infrared, iSwain et alJ (l2008T l used HST/NICMOS and 
showed evidence for absorption of atmospheric water and 
methane while Spitzer broadband transit p hotometry has 
given evidence for CO ([Desert et alj|2009al ). 

The current transmission spectrum of HD189733b ap- 
pears to be one which is dominated by condensed haze par- 
ticles in the optical, but dominated by the absorption of 
gaseous molecular species (e.g., H 2 0, CO, CH4) at near- 
infrared and infrared wavelengths. The transition wave- 
length between being dominated by haze or gaseous species 
should give indications of the particle size of the condensed 
haze, though it is not clear if all the current observations 
are completely compatible and can give a coherent picture 
of the planet's atmosphere. 

Here we present near-infrared HST photometry of 
HD 189733b during five transits, designed with the intended 
goal of measuring the signature of atmospheric water. This 
work is part of our ongoing efforts to charaterize the tran- 
sit spectra of HD18 9733 using space-based observatories 
(|Desert et al.ll200"9aU bT). We observed with two narrowband 
filters and designed the observing program to be unbiased 
from the typical orbit-to-orbit and intra-orbit systematic 
errors inherent in obs erving in spectroscopic modes with 



HST STIS, ACS (e. g.,lBrownet~aLl i2001: Pont et al. 2007) 



and NICMOS re.g.. lSwain et all 120081 : [Carter et al.ll2009h 
instruments. In Sect. 2 and Sect. 3, we detail these obser- 
vations, which include correcting for the presence of unoc- 
culted star spots and stellar limb darkening. In Sect. 4 we 
present the analysis of our transit light curves and discuss 
our results in Sect. 5 and Sect. 6. 
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Fig. 1. Aperture photometry from Visit #1 using the F187N 
narrowband filter. The data from all four spacecraft orbits of 
the HST visit are plotted back-to-back (there are ~45 minute 
data gaps between HST orbits due to Earth occultation). The 
transit occurs in orbit 3 between exposures 371 and 580 and is 
visible in the top plot. Plotted are the total counts from a 15.2 
pixel aperture (top), along with the counts from only the PSF 
central core (middle, 4.25 pixel aperture), and counts from only 
the PSF wings containing the 1 st and 2 nd airy rings (bottom, 
ring with inner-annulus of 4.25 pixels outer-annulus 15.2 pixels). 
The intra-orbit thermal focus variations change the shape of 
the PSF during an orbit, which can be seen by the large flux 
changes in counts of the PSF core and wings. However, with 
our narrowband filters the total counts are conserved in wide 
apertures, which exhibit no residual intra-orbit variations. A 
small intrapixel sensitivity related error is observed in the PSF 
core photometry which contaminates ~10% of the exposures 
(see text). 



2. Observations 

2.1. Hubble Space Telescope NICMOS narrowband 
photometry 

We observed HD 189733 during five transits using the 
Near Infrared Camera and Multi Object Spectrometer 
(NICMOS) aboard the Hubble Space Telescope during 
Cycle 16 (GO-11117). Each transit observation consisted 
of four consecutive spacecraft orbits, each roughly centered 
on a transit event. The visits occurred during: 15 April 2008 
UT 15:36 to 21:17 for visit #1; 3 May 2008 UT 08:51 to 
14:26 for visit #2; 18 May 2008 UT 22:44 to 19 May 2008 
UT 04:27 for visit #3; 25 May 2008 UT 14:33 to 20:15 for 
visit #4; and 11 Aug 2008 UT 06:24 to 12:04 for visit #5. 
For each transit event, we obtained images using the high 
resolution NIC1 camera with only a single narrowband fil- 
ter, using either the F166N filter centered at 1.6607 psm 
(AA = 0.0170/im) or the F187N filter centered at 1.8748 
fim (AA = 0.0191/im). By not switching between filters, 
the duty cycle of the instrument is increased as much less 
overhead time is spent between consecutive exposures. 

The NICMOS instrument has three 256x256 pixel 
HgCdTe actively cooled cameras with the NIC1 camera of- 
fering the highest spatial resolution, having a 11x11 arcsec 
field of view and 43 milliarsec sized pixels. The high resolu- 
tion of the camera combined with the use of narrowband fil- 
ters allowed us to observe the bright (V mag =7.67) star with 



short 2-4 second exposures well before saturating the detec- 
tor. We also adopted a modest defocus during our observa- 
tions, such that the photometry would be more insensitive 
to flat fielding errors and sub-pixel variations. The F187N 
filter images were taken in Multiaccum SPARS4 mode with 
Nsamp=3 yielding 4.00 second exposures while the F166N 
filter images were taken with Multiaccum STEP1 mode 
with Nsamp=4 yielding 1.99 second exposures. Both ob- 
serving modes have a 9 second readout overhead. In visit #1 
we obtained 793 F187N exposures, in visit #2 860 F166N 
exposures, in visit #3 418 F187N exposures, in visit #4 927 
F166N exposures, and in visit #5 751 F187N exposures. 
The data quality of visit #3 was compromised by a failure 
of a known 2-gyro issue with the Fixed Head Star Trackers. 
The failure caused a total loss of the visit #3 out-of-transit 
data. However, even without a long out-of-transit baseline 
flux, visit #3 still proved useful when analyzing the data 
in conjunction with the other two F187N visits. 

During each HST visit, we obtained 50 exposures in one 
filter during the beginning of the first orbit, before switch- 
ing to the opposite narrowband filter in which the rest of 
the visit, including the transit, was observed. This short set 
of opposite-band images was used to help monitor the stel- 
lar activity, such that the absolute flux and activity level 
of the star would be known at each of our two wavelengths 
for all five visits. 

We performed aperture photometry on the calibrated 
STScI pipeline reduced images. The pipeline includes cor- 
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rections for bias subtraction, dark current, detector non- 
linearities, and applies a flat field calibration. The aperture 
location for each image was determined using a two dimen- 
sional Gaussian fit of the point spread function (PSF) for 
each image. We found apertures with a radii of 15.2 pix- 
els for the F187N filter and 14.7 pixels for the F166N filter 
minimized the standard deviation of the out-of-transit light 
curves. The background in our short exposures was found 
to be negligible, typically accounting for only ~50 total 
counts per image. 

Our choice of a wide aperture minimizes the effect 
of intra-orbit thermal focus variations (see Fig. 1). This 
well known effect (see STScI Instrument Science Reports 
ACS2008-0i] and NICMOS ISR2007-0033) changes the 
shape of the PSF during an orbit and is due to small 
changes in the focus of the telescope, which are ascribed 
to temperature variations. With a sufficiently wide aper- 
ture, covering at least the central peak and first and second 
airy rings, no residual intra-orbit variations are observed. 
In four of our five orbits, the photometric time series ex- 
hibit sudden drops in flux, confined to the central core of 
the PSF and effects 10-15% of the images. We attribute 
this phenomenon to a likely intrapixel sensitivity variation, 
as it is confined to pixels in the central PSF core which 
have the largest change in flux across a pixel and the ef- 
fect is smaller for higher levels of defocus. Images which 
show this intrapixel sensitivity effect were easily flagged by 
comparing the photometric time series of the photometry 
of the central PSF core to the photometry of the outer-airy 
rings, which are insensitive to the effect (see Fig. 1). We 
adopted two strategics to deal with these flagged images, 
choosing to either eliminate them from the remaining anal- 
ysis or interpolate the affected PSF core photometry using 
neighboring images in the time series. Both methods pro- 
duced equivalent results in the final transit light curve fits. 
We note that future HST/NICMOS programs can eliminate 
this small problem completely by adopting a larger defocus, 
as is the case here in visit #5 which does not exhibit this 
effect. 



2.2. Monitoring the stellar activity 

We monitored the stellar activity of HD 189733 using both 
ground-based data and the absolute flux level from the 
NICMOS instrument itself. The ground-based coverage 
was provided by the T10 0.8 m Automated Photoelectric 
Telescope (APT) at Fairborn Observatory in southern 
Arizona. This ongoing observing campaig n of HD189733 
began in October 2005 and is detailed in iHenrv fe~W mn 
(|2008T ). The APT uses two photomultiplier tubes to simul- 
taneously gather Stromgren b and y photometry. The APT 
data spans four of our five HST visits (see Fig. 2), with only 
the fifth visit not covered as the telescope is shut down dur- 
ing the Arizona monsoon season. This dataset also covers 
the HD189733 stellar activi ty during the epoch s of the HST 
trans mission spectra from iPont et alj (l2007f) and Swain 
et al. ( | 2008r) as well as the Spitzer transit photometry of 
IDesert et ap (l2009af). IDesert et al.l (l2009bf). and Knutson 
et al. (2009) making it an invaluable resource when com- 
paring the stellar activity level from epoch to epoch. 




http: / / www.stsci.edu /hst /acs/documents / isrs / 
2 http:/ /www. stsci.edu/hst/nicmos/documents/isrs/ 



Fig. 2. Simultaneous ground-based pho tometric monitoring o f 
HD189733 (blue crosses) continued from lHenrv fc Winnl {2008), 
along with a fitted sinusoid (blue). The NICMOS flux levels 
for all five visits are also plotted for both the F166N filter (red 
squares) and the F187N filter (red dots). Assuming 1000 K cooler 
star spots with ~ 2% coverage, the amplitude of the near-IR 
spot modulation is expected to be about half that observed in 
the optical (red sinusoid). The NICMOS flux levels are seen 
to modulate in phase with the optical ground-based data at 
approximately the expected near-IR spot amplitude. 



The near-IR flux levels for HD189733, during our five 
visits, are based on the baseline flux level from each transit 
light curve as well as the first 50 exposures obtained in the 
filter opposite to that of the transit (ie. F166N or F187N). 
While these first 50 exposures were obtained during the 
first orbit of a visit, these orbits did not show any signifi- 
cant signs of systematic errors. Our use of the first orbit is 
in contrast to nearly all previous similar HST studies (in- 
cluding the STIS and NICMOS instruments), which have 
had to disregarded the first orbit of an HST visit. During 
the first orbit, the telescope thermally relaxes in its new 
pointing position causing significant flux variations. The 
near-IR stellar flux levels are shown in Fig. 2, with an arbi- 
trary zero-magnitude flux, which clearly modulate in phase 
with the optical measurements. 

3. Stellar limb darkening 

3.1. Three parameter non-linear limb-darkening law 

For solar-type stars at near-infrared and infrared wave- 
lengths, the strength of stellar limb darkening is weaker 
compared to optical wavelengths. However, the intensity 
distribution is increas ingly non-linear a t these longer wave- 
lengths (see Fig. 5 of ISing et al.ll2008a! ). which can require 
adopting non-linear limb-darkening laws when fitting high 
S/N transit light curves. Many studies have chosen to use a 
quadratic or four parameter non-linear limb- darkening law 
to describe the stellar intensity distribution ( Claret 2000), 
with the coefficients determined f r om ste llar atmospheric 
models. As noted by ISing et ail (|2008aD for the widely 
used Kurutz ID ATLAS stellar models, the largest dif- 
ferences between existing limb-darkening data and the ID 
stellar models is at the very limb, where ATLAS models 
predict a dramatic increase in the strength of limb dark- 
ening. For the sun, the ATLAS models over-predict the 
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Table 1. Three parameter non- linear limb darkening coef- 
ficients 



coefficient 


Ff66N filter 


F187N filter 


C2 


2.1483 


1.9508 


C3 


-2.7763 


-2.4507 


C 4 


1.1265 


0.9770 



strength of limb darkening by >20% at fi = cos{6) values 
below 0.05, though they per form well otherwise. Both solar 
data (iNeckel & Labslll994ri and 3D stellar models ( Bigot 
. 1200 " 



et al. 2006) show the intensity distribution at the limb to 
vary smoothly to ji = 0, with no dramatic or sudden in- 
creases in limb-darkening strength. To mitigate this limb 
effect, we choose to fit for the limb-darkening coefficients 
from ATLAS models using only values of fi > 0.05. We also 
choose to fit a three parameter non-linear limb-darkening 
law, 



1 - c 2 (l - /x) - c 3 (l - M 3/2 ) - c 4 (l - M 2 ) 



(1) 



m 
i{±) 

as the /j, 1 / 2 term from the four parameter non-linear law 
mainly affects the intensity distribution at small /i values 
and is not needed when the intensity at the limb is de- 
sired to vary linearly at small [i values. Compared to the 
quadratic law, the added /i 3 / 2 term provides the flexibil- 
ity needed to more accurately reproduce the stellar model 
atmospheric intensity distribution at near-infrared and in- 
frared wavelengths. 




2.0 3.0 
Wavelength fj.m 



Fig. 3. The impact parameter (top) and system scale parameter 
(bottom) vs. wavelength for HD189733 as measured by Desert 
et al. ll2009d).lDesert et al. l ll2009bf).lKnutson et all (12009D, Pont 
et al. (|2007l ). I Winn et al.1 l|2007h . and this work. Excluding two 
deviant Spitzer valu es from Desert et al. (2009; see also Paper 11 
iDesert et al. 2009b for further details), there is a good general 
agreement between the various HST, Spitzer, and ground-based 
results with a mean 6 value of 0.6636±0.0031 and mean a/R* 
value of 8.898±0.027. These mean values are over-plotted (hor- 
izontal line) along with their la errors (dashed vertical lines) 
and adopted for our final transit light curve fits seen in Figures 
5 and 6. 



3.2. HD189733 limb darkening 

The limb-darkening coefficients for the three parameter 
non-linear law were computed using a Kurucz ATLAS stel- 
lar mode0 with T e// =5000 K, log g=4.5, and [Fe/H]=0.0 
in conjunction with the transmission through our two nar- 
rowband filters. For reasons described above, we fit for 
the limb-darkening coefficients using the calculated inten- 
sities between /i=0.05 and (x=l, providing the coefficients 
in Table 1. 



4. Reduction and analysis 

4.1. Correcting for non-occulted stellar spots 

A transit light curve is affected by the presence of star 
spots, both when a spot is occulted by the planet and by 
the presence of non-o cculted spots duri ng the epoch of the 
transit observations (|Pont et al.ll2008l ). We will focus on 
non-occulted spots here, as there is no evidence for occulted 
spots in any of our observed transit light curves. 

Non-occulted stellar spots affect the shape of a tran- 
sit light curve by making the stellar surface dimmer than 
a spot-free stellar disk. During transit, the planet hides a 
larger fraction of the total overall flux, thereby increasing 
the apparent size of the planet. A wavelength dependance 
is also introduced, as the stellar flux lost by the spots at 
a given wavelength will depend on the blackbody temper- 
ature difference between the stellar surface and the spots. 
The stellar spots also create a quasi-periodic photometric 
variability, as the spots rotate into and out of view at the 
stellar rotation period. 



http:/ /kurucz. harvard.edu/ 



To correct for the presence of non-occulted stellar spots 
in this program, we used the ongoing ground-ba sed photo- 
metric monitoring data of iHenrv fe: Winn] (|2008t ) as well as 
the calibrated flux level of the N1CMOS observations, mea- 
suring the flux variation due to stellar spot modulation. 
HD189733 is an active star, varying in the optical at the 2- 
3% level with spot-coverage at any given epoch covering ~1- 
2% of the stellar surface. The NICMOS observations vary 
by ~1% (see Fig. 2), which is consistent with the expected 
near-IR variation due to stellar spots as it occurs both in 
phase with the ground-based photometric monitoring and 
near the amplitude expected w hen assuming ~40 00 K spot 
temperatures as measured by iPont et al.l ( 20081 ) and the 
observed optical spot coverage. Furthermore, for each visit 
both the F166N and F187N filters showed similar near- 
IR stellar flux variations, as would be expected for stellar 
spots. The baseline flux of each visit was corrected in a 
differential manner, adding the necessary flux to each visit 
such that all the observations are compared at the same 
minimum spot activity baseline flux level. HST visits ^=1 
and #5 occurred at a spot activity minimum, and were 
chosen here as the flux reference for the remaining visits. 
The corresponding effect on the final determined radius is 
somewhat small but these corrections are important as ex- 
oplanetary atmospheric signatures are typically also small. 
An estimated ~0.0015 uncertainty on the absolute flux level 
of each visit translates into a spot-correction uncertainty of 
~0.0001 Rpz/R* and ~0.25tr on the final values of R p ;/R* 
for the two wavelengths when combining the different visits. 
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4.2. Transit light curve fits 

We modeled the transit light curve w ith the theoretical 
transit models of iMandel fc Ag ol (2002). We choose to fix 
the planetary orbital phase using the ephcmcris of Knutson 
et al. ((20091) (P = 2.21857578 ± 0.00000080 days; Tc = 
2454399.23990 ± 0.00017 HJD), fitting for the planet-to- 
star radius contrast R p ;/R*, the inclination i, and the mean 
stellar density , which is proportional to the system scale 
(a/R+) cubed. No significa nt deviations from the Knutson 
et al. (|2009l ) ephemeris were found, with the three F187N 
visits determined to be within 2 seconds of the predicted 
central transit time. To account for the effects of limb- 
darkening on the transit light curve, we adopted the three 
parameter limb-darkening law and let the linear coefficient, 
c 2 , free in the fit, while C3 and C4 were fixed to the best- 
fit model values listed in Table 1. Choosing to fit only one 
limb-darkening parameter ensures that the fits are both not 
significantly biased by adopting a stellar atmospheric model 
and do not suffer from degeneracies between fitting multiple 
limb-darkening parameters. We also allowed the flux level 
of each visit to vary in time linearly, described by two fit 
parameters. The linear trend accounted for any possible de- 
tector drifts or stellar activity related flux variations during 
the observations and we found that visits # 4 & 5 exhib- 
ited a small slope. The best fit parameters were determined 
simultaneously with a Levenberg-Marquardt least-squares 
algorithm using the unbinned data (see Table 2), with the 
error of each datapoint set to the standard deviation of the 
out-of-transit residuals of a given visit. 

The fit values for the linear limb-darkening coefficients 
are consistent within la of the model values (see Table 2). 
We also performed fits letting either C3 or C4 free, with C2 
fixed, finding similar results indicating that the Atlas mod- 
els are performing well. As a check of our choice adopt- 
ing the three parameter law, we also fit the four-parameter 
limb-darkening law fixing C2, C3, and C4 to their model val- 
ues and letting the c\ coefficient (proportional to p 1 ^ 2 ) free. 
At both 1.66 and 1.87/im, we found c\ to be consistent with 
zero. In addition, we tried fits with the quadratic limb- 
darkening law fi tting both coefficients in a minimally cor- 
related fashion (|Pall 12008). finding equivalent transit pa- 
rameters, though the three-parameter law still produced 
the best x 2 fits. 

Given the accumulating high precision transit data on 
HD189733 with both HST and Spitzer, several orbital sys- 
tem parameters (such as i, b, and a/R*) can be determined 
more precisely than with just these HST/NICMOS transit 
observations alone. Transit-fit values for R p ;/R^ are cor- 
related to these system parameters, which can adversely 
affect a proper comparison of the planetary radii across 
different wavelengths, when searching for atmospheric sig- 
natures, if significantly deviant values are used. Therefore, 
it is advantageous to fit for R p ;/R* using the best system 
parameters available, and also necessary that the param- 
eters are consistent across different studies when compar- 
ing the planetary radii between multiple observations. In 
Figure 3, we plot the measured orb it al system pa r ameter s 
of HD189733 f romlPont et alj (120071) . Desert et all (l2009af). 
iDesert et all (|2009bf ). iKnutson et all (j2009l ). IWinn et all 
(|2007f ). and this work using the individual NICMOS fits. 
The results for i, a/R*, b, and nearly all agree at the 
~ lcr level. A nota ble exception i s the 3 .6/nro and 5.8fim 
Spitzer results from lDesert et ah! (|2009af ). which show sig- 



0.006 
0.004 
0.002 




Fig. 4. NICMOS filter transit photometry of all three visits 
using the 1.8748^m narrowband filter. The best-fit transit 
light curve model is also fit (grey) along with the observed - 
calculated residuals of the fit (bottom panel). The light curve 
has a S/N of 1000 per-point and an equivalent precision of 
4.0xl0~ 4 in an 80 second bin. 
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Fig. 5. NICMOS filter transit photometry of both visits using 
the 1.6607/xm narrowband filter. The best-fit transit light curve 
model is also fit (grey) along with the observed - calculated 
residuals of the fit (bottom panel). The light curve has a S/N of 
600 per-point and an equivalent precision of 6.3 xlO -4 in an 80 
second bin. 



nificantly deviant values in both b and a/R*. These deviant 
values may have been the result of occulting stellar spots 
during ingress or egress and were disregarded in this study 
(see Desert et al. 2009b for further details). From the HST, 
Spitzer, and ground-based studies, we derived mean orbital 
system parameter values calculating a weighted mean using 
the different measurements. The determined mean values 
and their associated error are listed in Table 2. These val- 
ues have an improved determination of i, 6, and a/R* by 
up to a factor of seven compared to our individual NICMOS 
fits and were adopted for our final transit light curve fits 
(labeled Joint HST & Spitzer in Table 2 and plotted in Fig. 
4 and 5) . We also choose to adopt the model limb-darkening 
values from Table 1 for the three coefficients, minimizing 
the number of free fit parameters. 
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Table 2. System parameters of HD189733b 



Pa T-a TTIP^PT 


Individual NICMOS Fits 


Tnint HST h Snit7Pr T 




r loolN niter 


r 18 / IN niter 


Flfir-lN filtpr F1&7N filfpr 

1UU1N lllLtJI r 1CW1N lllLcI 


central wavelength, A (/im) 


f.6607 


f.8748 




wavelength range, AX(pm) 


0.0f70 


0.0191 




planet-to-star radius contrast, R p ;/R* 


0.f5464±0.0005f 


0.15496±0.00028 


0.15498±0.00035 0.15517±0.00019 


inclination, i (deg) 


85.75±0.f3 


85.715±0.076 


85.723±0.024 f 


system scale, a/R* 


8.87f±0.095 


8.854±0.062 


8.898±0.027 f 


impact parameter, 6=ocosj/R* 


0.657±0.02f 


0.661±0.013 




stellar density, (g cm -3 ) 


2.684±0.086 


2.669±0.056 


2.709±0.025 t 


linear limb-darkening coefficient, ci 


2.f95±0.058 


1.972±0.029 


2.1483 1.9508 



^System parameter values and 1-<t error derived from the weighted mean using: Pont et. al. (2007) HST/ACS, Desert et 
al. (2009a,b) Spitzer/IRAC, Knutson et al. (2009) Spitzer/IRAC, and the individual fits of this work HST/NICMOS. The mean 
value was adopted for the fits. 



4.3. Red noise estimation 

The RMS of the residuals between the best-fitting Joint 
HST & Spitzer models and the data are 9.98 xlCP 4 for the 
F187N filter and 1.67xl(T 3 for the F166N filter. These val- 
ues are approximately 1.5 times the calculated pipeline er- 
ror, which is dominated by photon noise. This is a modest 
improvement over other similar NIC MOS transit studies 
((Carter et al.ll2009l ; IPont et al.ll2009l ). which have gener- 
ally found 2-3 times the photon noise limit. We checked 
for the pres ence of systemat ic errors correlated in time 
("red noise" . IPont et al.ll2006D using two methods. We first 
checked that the binned residuals followed a iV -1 / 2 relation 
when binning in time by N points. The presence of red noise 
causes the variance to follow a a 2 = cr^/N + <r 2 relation, 
where a w is the uncorrelated white noise component while 
ay characterizes the red noise. We found no significant ev- 
idence for red noise in either the F187N data or F166N 
data, when binning on timescales up to the ingress and 
egress duration. We also used the "prayer-bead" method to 
check for red n oise, which consists o f a residual permuta- 
tion algorithm (|Moutou et al.ll2004D . In this method, the 
residuals of the best-fit model are shifted with wraparound 
over the light curve and a new fit is performed. The struc- 
ture of any time correlated noise is preserved and its effect 
on the uncertainty of the transit light curve parameters are 
revealed in the distribution of fit parameters. The prayer- 
bead method also showed no significant signs of red noise, 
as the distribution of each of the fit parameters at both 
wavelengths were compatible with the uncertainties found 
from standard x 2 statistics. With no detectable red noise, 
our S/N levels are comparable to other high S/N transit 
light curves, with precisions of 4.0xl0 -4 per 80 second bin 
for F187N and 6.3xlCT 4 per 80 second bin for F166N. 

5. Discussion 

5.1. Expected H2O atmospheric signature 

We have presented high signal-to-noise near-infrared tran- 
sit narrowband photometry of HD189733 at 1.66 and 1.87 
/im, with the intent of searching for a transit signature of 
atmospheric water. If atmospheric H2O were the dominant 
absorbing species at both of these wavelengths, we could 
expect significant radii differences as H2O has a strong 
absorption bandhead over the probed wavelength range. 
The expected difference in radii can be estimated using 
the difference between opacities at 1.66/irra and 1.87/im. 



We calculate the effective absorption cross sections per 
molecule to be cti. 66 =8.54x 10~ 24 and cti. 8 7=1-79x10~ 21 
(cm 2 mol ecule" 1 ), resp e ctivel y, using the high temperature 
line list of lBarber et al.l (|2006t l , a Voigt broadening line pro- 
file, a temperature of 1200 K, and the response through our 
two narrowband filters. The variation in apparent planetary 
radius then follo ws from Eg. 2 of ( Lecavelier des Etangs 
et al. I2008al ) giving. 

AiZj, = iTm — , (2) 
ci.66 

where H is the atmospheric scale height. H is given by H = 
kT / p,g, where u is the mean mass of the atmospheric par- 
ticles taken to be 2.3 times the proton mass, T the temper- 
ature, and g the surface gravity. Assuming i?* = 0.766 Rq, 
R p i = 1.155 Rjup, M p i = 1.144 Mj up , and a temperature 
of 1200 K, we find H/R^—Q. 00038 which gives an expected 
radius difference of 0.0020 R p/ /R*. Our HST NICMOS ob- 
servations show only a radius contrast of 0.00019±0. 00040 
Rpz/R* between the two wavelengths, a factor of 10 lower 
and 4.6(7 away from the expected value. If H2O was the 
dominant opacity source in HD189733b's atmosphere at 
these altitudes, it would have been easily detected by our 
observations. 

5.2. Comparisons with previous results 

Swain et al. ((20081 ) used HST NICMOS to produce a 
near-IR low resolution grism transmission spectrum of 
HD189733b, which appears to show two strong atmospheric 
features at H2O and CH4. This grism spectrum encom- 
passes our wavelength range (see Fig. 6) and appears to 
show a strong H2O absorption molecular bandhead feature 
betw een 1.5 and 2.0 um. Between 1.66 and 1.87 am. Swain 
et al. ((200^) found a radii contrast of 0.00233±0. 00035 
Rpj/R*. This radii contrast is significantly larger than the 
one we observe, with our results failing to reproduce the 
Swain et al. value at a 5cr confidence level. 

Given th e large disagreement between our results and 
iSwain et all (2008), we find two scenarios that could pos- 
sibly explain the conflicting results. Either the planet's 
transmission spectrum is variable, or residual systematic 
error s still plague the ed ges of the Swain et al. spectrum. 
Both lSwain et al.l (|2008l ) and our results correct for stellar 
starspot activity, ruling out stellar variability as the under- 
lying cause. In addition, the discrepancy between the two 
results is too large to attribute to random statistical fluctu- 
ations. Variability in the planet's transmission spectrum by 
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0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Wavelength (/xm) 

Fig. 6. Measured planetary radii for HD189733b at optical and near-infrared wavelengths. Plotted are the results from our 
NICMOS narrowband photometry (black dots), along with the NICMOS grism spectrum (red squares) from Swain et al. (2008), 
and ACS grism spectrum (green dots) from Pont et al. (2008). The l-cr error bars on the fit radii are indicated (y-axis error bars), 
along with the wavelength range of each observation (x-axis error bars, grey vertical bars). Also plotted (purple) is the prediction 
by Rayleigh scattering due to haze from Lecavelier et al. (2008a), projected here into the near-infrared along with the 1-a error 
on the predicted slope (purple, dashed lines). The NICMOS spectrum from Swain et al. (2008) is quoted as being uncertain in 
its absolute flux level by ±2xl0" 4 (±0.00064 R pi /R*), an offset of -0.00042 R p ;/R* was applied here for comparison reasons such 
that the values at 1.87 /im match. Our 1.66/j.m results are in disagreement with the both the Swain et al. spectra and the expected 
H2O atmospheric signature (blue dots), but are in excellent agreement with the predicted planetary radii values from atmospheric 
haze. 



dynamic weather processes would seem a plausible scenario, 
especially given the seemingly variable emission spectrum 
results from Spitzer (jGrillmair et al.ll2007l . [20081 ) . However, 
we find the transmission variability scenario difficult to sub- 
stantiate, particularly with the large systematic errors in- 
herent with NICMOS grism observations, which can com- 
promise precision photometric work. In order for the vari- 
ability scenario to be viable, the optically thick absorbing 
haze would have to clear out of the atmosphere, over the en- 
tire limb, down to altitudes at least 6.7 scale heights lower 
(~1200 km) to reveal the full extent of the H2O feature. 
The haze would also have to clear for only the single epoch 
of the Swain et al. (2008) HST visit, but not during the 
seven epochs of the HST visits between Pont et al. (2008) 
and this work. Likely, the variability scenario is also de- 
pendent upon the haze composition, and the transmission 
spectrum would be difficult to vary with non-condensate 
Rayleigh sc attering candidates such as H? ( Lecavelier des 
Etangs et al. l2008al ). We note that in our HST visits span- 
ning five months, no significant visit-to- visit planet radius 
differences were observed at either wavelength, after cor- 
rection for non-occulted stellar spots. At 1.66 fim we find 
radii of 0.1545±0.0004 and 0.1551±0.0006 R p/ /R* for vis- 



its #2 and 4, respectively, while at 1.88 /im we find radii 
of 0.1549±0.0003 and 0.1558±0.0004 R p; /R* for visits #1 
and 5, respectively. 

HST NICMOS slit-less grisms have a number of instru- 
ment related systematic effects that appear in precision 
photometry. These effects include detector wavelength de- 
pendancies such as the flat field, sensitivity, PSF, and reso- 
lution which interact in a complex manner with telescope's 
PSF stability and pointing accuracy. With these effects, 
typical transit programs with HST NICMOS grisms tend 
to exhibit large intra-orbit flux variations as the telescope 
thermally relaxes from day-to-night and orbit-to-orbit flux 
variations largely dependent on the grating fi lter wheel po- 
sition and spectral trace (see Fig. 1 & 2 of ICarter et al.l 
120091 ) . In addition, the low resolution Nic 3 camera used 
for grism observations is highly under-sampled leading to 
large intra-pixel sensitivities, though these can be largely 
averaged out by defocusing the detector. The importance of 
systematic errors resulting from the grating filter wheel po- 
sition on HD18 9733b's spectrum can be seen in the emission 
spectrum from ISwain et aTl ([2009). This spectrum shows 
nonphysical results (with negative planet fluxes) at the 
blue-edge of the spectrum. At the edges, the spectral trace 
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has the largest orbit-to-orbit deviations, causing greater 
systematic errors. The HD189733b NICMOS transmission 
and emission spectra were both observed and reduced in 
much the same manner by Swain et al., thus the edges 
of the transmission spectra should also be suspect of un- 
derestimated uncertainties and greater residual systematic 
errors. 

In a separate NICMOS spectroscopic study of GJ436b 
by Pont et al. (2008), variations at the 4 x 10 -4 flux level in 
the transmission spectrum were also found, similar in am- 
plitude to the Swain et al. spectral features of HD189733b. 
However, given the level of measured systematic errors, 
Pont et al. (2008) chose not to interpret the GJ436b vari- 
ations as due to real planetary transmission features, but 
rather as systematic errors. These errors exhibited a wave- 
length dependence that was time dependent, with spec- 
tral features slowly shifting in wavelength with time, as 
expected from the HST thermal focus variation systematic 
effect. While Pont et al. (2008) and Swain et al. (2008) used 
different NICMOS grisms, both corrected for systematic 
errors in a similar manner adopting a multilinear decor- 
relation against external variables including the position, 
rotation, and width of the spectral trace. Thus, similar sys- 
tematic errors between the two studies may be expected. 
A systematic error of 4 x 10~ 4 in flux would correspond 
to an error of ^0.0013 Rpj/R* at a given wavelength in 
HD189733b and ~0.0018 R p ;/R* when comparing the dif- 
ference between two wavelengths, possibly accounting for 
the features in the Swain et al. (2008) spectrum. 

With NICMOS grisms, wavelength shifts and a variable 
PSF with time can change the flux per pixel with time, 
which becomes compounded by pixel-to-pixel and sub- 
pixel response variations which are not accurately known. 
Observing with narrowband photometric filters mitigates 
these wavelength dependent and PSF effects, eliminating 
the need to correct for the systematic errors shown by the 
grism spectra (a main motivation behind our original selec- 
tion of narrowband filters). Given our robust (decorrelation 
free) photometry and the suspect edges of the t ransmission 
spectr a, the differences between our results and lSwain et al.l 
( 2008) would clearly seem to favor grism-related systematic 
errors as the principal cause. 

5.3. Atmospheric haze in the near-infrared 

Although our results on the planet-to-star radius ratios 
measured at 1. 66 /im and 1 . 87 are significantly dif- 
ferent from the ISwain et alj (|2008l ) measurements at the 
same wavelengths, our results are still consistent at both 
wavelengths with the extrapolation of measurements at op- 
tical wavelengths assuming R ayleigh scat t ering . In effect, 
the accurate measurements of iPont et al.l ()2008f ) showed a 
significant decrease of the planet radius toward longer wave- 
lengths t hat was interpreted as the signature of absorption 
by haze. iLecavelier des Etangs et alT (p008a') showed that 
the radius variation from 0.55 to 1.05 /im is a signature of 
the Rayleigh scattering with a decrease of the cross-section 
following a A power law. This law can be extrapolated to 
1.87 /im. We find that both our measurements at 1.66 /im 
and 1.87 /im are in excellent agreement with the Rayleigh 
scattering predicti on (Fig. 6V As shown in Lecavelier des 
Etangs et al. (2008al), the slope of the radius as a function 
of the wavelength in the Rayleigh regime is characteristic of 
the atmospheric temperature. The slope for extreme tem- 



0.1530 




2.0 2.2 
Wavelength (/im) 



Fig. 7. Similar to Fig. 6, except the Swain et al. (2008) data 
are plotted at the original unshifted level and the points at the 
ends of the spectrum (likely affected by systematic errors) are 
removed. Both lcr error-bars (thick) and 2cr error-bars (thin) on 
the Swain et al. spectra are also shown to illustrate that it is 
consistent at that level with Rayleigh scattering by condensate 
haze, below 2 fim. 



peratures of 1190 and 1490 K are plotted in Fig. 6 (dotted 
purple lines). We can see that the smaller slope correspond- 
ing to the lowest temperatures is slightly favored, in par- 
ticular by the 1.87 /im measurement. 

We performed a glob al fit to the presen t new measure- 
ments together with the IPont et all (|2008f ) published val- 
ues. With the hypothesis of pure Rayleigh scattering, the 
data are well fitted with a reduced \ 2 of 1.1 despite the 
extended width of the spectral domain spread over a factor 
4 in wavelength. We find that the temperature given by the 
Rayleigh power law is 1280±110K. This temperature agrees 
with Spitzer 8 and 24 fim phase curve measurements, which 



also find similar temper atures of (T„ 



-T mm )/2 = 1134 



and 1102 K respectively (jKnutson et al.|[2~009T ). Relaxing the 
Rayleigh hypothesis and using Mie scattering from 0.55 /im 
to 1.87 /im, we obtain new con straints on the particle size 
and th eir optical properties (see lLecavelier des Etangs et alj 
l2008af ) , with the fit to the data slighly improved. Mie scat- 
tering implies a larger possible temperature range of 1170 
to 1870 K (includes lcr error) because of a less steep extinc- 
tion variation with wavelength in the near infrared. The 
proposed scenario with MgSi03 still holds, and we find that 
the particle maximum size must be in the range 0.009 /im 
to 0.086 /im. Recall that because the Rayleigh scattering 
cross-section is proportional to the particle size, a par t, to 
the power of 6 (cr oc ctp art ), the size range given above is the 
size of the largest particles in the size distribution. 

Finally, the interpretation of the present transit absorp- 
tion signatures in terms of Rayleigh scattering at A 5s 2 /im 
is consistent with the measurements of emission spectra 
which do not show signatures of haze absorption. This is 
explained by the fact that the slant transit geometry probes 
atmospheres at high altitude and is more sensitive to low 
absorption cross-sections, while the normal viewing geome- 
try of emission spectra probes deeper in the atmosphere 
()Fortnevl 12005). Though more uncertain in the absolute 
planet-to-star radius contrast, the lSwam" et all (2008) trans- 
mission spectrum can also be considered as consistent with 
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the haze interpretation if the extreme parts of their absorp- 
tion spectra are disregarded (Fig. 7). 

6. Conclusions 

Narrow-band NICMOS photometry can provide a robust 
method for obtaining precision near-infrared photometry. 
With transiting planets, this can lead to secure near- 
infrared planetary radii where important molecular and 
condensate atmospheric planetary signatures can be found. 
With two carefully selected narrow-band filters, we are able 
to provide stringent constraints on the presence of H2O 
absorption features. We find that upper-atmospheric haze 
provides sufficient opacity in slant transit geometry to ob- 
scure the near-IR H2O molecular signatures below 2/jm. 
While absorption due to H 2 was reported by I Swain et al.1 
(12008ft . we rule out any such feature at 5er and find that 
their observed feature is likely an artifact of residual sys- 
tematic errors comprom ising the spectral e dges. The plan- 
etary radii values from iPont et alj (|2008l ) and this work 
indicate that Rayleigh scattering dominates the broadband 
transmission spectrum from at least 0.5 to 2 [im, with 
sub-micron MgSiOs haze particles a probable candidate 
(|Lecavelier des Etangs et all [200 8a). The transition wave- 
length between the transmission spectra being dominated 
by haze particles and molecular features would seem to 
be between 2 and 3 fim, as the IRAC Spitzer photome- 
try between 3 and 8 microns all show planet radii in ex- 
cess of those predicted by Rayleigh scattering ([Desert et al.l 
I2009al lbl). To be clear, while' our observations rule out de- 
tectable quantities of H2O in the near-IR at these high alti- 
tudes, the presence of methane in the Swain et al. transmis- 
sion spectrum is still possible, though it needs verification, 
as the 2.2 to 2.4 \im region appears to show an excess ab- 
sorption over the Rayleigh prediction. 

Spectroscopic characterization of transiting planets 
with JWST NIRCam and NIRSpec are antici pated to make 
signifi cant contributions to exoplanet science fD eming et al.l 
2009). These instruments have similar detectors and observ- 
ing modes as HST NICMOS, with NIRCam also containing 
similar narrow-band filters to those used here. While the 
systematic errors observed with HST and Spitzer should 
be greatly reduced in JWST, the correspondingly greater 
precisions that will be attempted, as super-Earths are ob- 
served, ensure that at some level these systematic errors 
will likely remain. Strategic use of narrowband filters can 
provide robust measurements with which to compare atmo- 
spheric features seen in transmission or emission spectra, 
checking both observing methods at key wavelengths. As 
no upcoming instrument will soon be able to compete with 
JWST, such independent methods may be the only way 
to confidently verify groundbreaking results in super-Earth 
spectra requiring significant systematic error corrections. 
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